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We present an attempt to build up a new two-body effective potential for hydrogen fluoride, fitted 
to theoretical and experimental data relevant not only to the gas and liquid phases, but also to the 
crystal. The model is simple enough to be used in Molecular Dynamics and Monte Carlo simulations. 
The potential consists of: a) an infra-molecular contribution, allowing for variations of the molecular 
length, plus b) an mier-molecular part, with three charged sites on each monomer and a Buckingham 
"exp— 6" interaction between fluorines. The model is able to reproduce a significant number of 
observables on the monomer, dimer, hexamer, solid and liquid forms of HF. The shortcomings of 
the model are pointed out and possible improvements are finally discussed. 
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I. INTRODUCTION 



Hydrogen fluoride (HF) is one of the simplest molecules capable of forming hydrogen bonds. Despite the simplicity 
of the compound, the theoretical description of its behavior, especially in the liquid phase, is still far from being 
fully satisfactory. Much relevant theoretical work has concentrated on the problem of determining a potential model 
suitable for computer simulations. For this purpose a model is needed that reproduces correctly the main features of 
the real interaction potential and which is simple enough to be computed efficiently. 



In the last two decades several emgjigcal potentials have been developed for Mjjl|jcular Dynamics (MD) or Monte 



iy) ' Carlo (MC) simulations of liquid HFBa Of these models, only the three-site one are able to reproduce correctly 
the dipole and quadrupole moments of the monomer. These models represent the charge distribution of each monomer 
by fractional charges placed at three sites on the molecular axis, two on the F and H nuclei and the third one at an apt 
propriate position X along the F-H bond. A first three-site model (called HF3) was developed by Klein andMcDonaldtJ 
by fitting ab initio results for the potential energy surface of the (HF)2 dimer. Cournoyer and Jorgenseria proposed a 
second three-site model (called HFC), with a simplified non-Coulombic part consisting of a single Lennard- Jones inter- 
action between the fluorines. The parameters of the model were fitted directly to experimental thermodynamic data 
for the liquid, n while simultaneously providing a reasonable equilibrium geometry of the dimer. Recently, Jedlovszky 
and Vallauritrel presented two further models, hereafter referred to as HF-JV1 and HF-JV2, respectively. HF-JVI is a 
variant of HFC, with charges reproducing both the monomer dipole and quadrupole and with an accurate treatment 
of the long range part of the Coulombic interactions, neglected in the original work on HFC.B The HF-JV2 model 
includes molecular polarizability, by adding an induced point dipole moment at each F site, while keeping the charge 
distribution of HF-JVI unchanged. The scalar polarizability of the molecules was set equal to its experimental value, 
while the two parameters of the Lennard-Jones interaction between F atoms were fitted to the experimental values of 
the liquid density and internal energy. 

Unfortunately, none of the available models is able to reproduce, in a fully satisfactory way, both the thermodynamics 
and the structure of liquid and gaseous HF. In the search for a new potential, the essential ingredients can be identified 
by reviewing the known properties of the gas, liquid and solid phases of HF. The simplest associated form of HF is the 
gas phase dimer, (HF)2, whose stipjjcture and rovibrational spectrum were first characterized by the microwave spectra 
of Dyke, Howard and KlemperermL2l The equilibrium configuration of the isolated dimer F-H' • • -F-H^ i*. planar but 
bent, as a consequence of a competition between dipolar and quadrupolar electrostatic interactionst2nl3 (whenever 
ambiguity is possible, H' denotes the hydrogen atom involved in the hydrogen bond, while H" is the other one). The 
atoms F-H' • • -F form a nearly linear arrangement, with H' placed slightly off the FF axis, and with a distance between 
the hydrogen-bonded fluorines tff ~ 2.7 A. The second hyfthftgen atom forms an angle FFH" « 115° and the H-F 
bonds are slightly stretched with respect to the monomer.il 



1 



The (HF)2 bent dimer appears to be the basic structural motif of-all the best known associated forms of HF, namely 
the gas phase (HF)@ cyclic hexamer Ji3 the low temperature crystaO and even the liquidO Pairs of adjacent HF units 
in the hexamer exhibit a bent arrangement similar to that of the dimer. The same structure is found in crystalline 
deuterium fluoride (DF) which is made up of infinite zig-zag chains of DF units£3 Interchain F-F distances (« 3.2 A) 
are much larger than intrachain F-F distances (~ 2.5 A), an indication that interchain interactions are weaker than 
the interactions, between adjacent HF units within the chains. The very small entropy difference between,-the solid 
and the liquicfca suggests that the liquid is largely associated. Best fit analysis of dielectric constant datallj and of 
Raman and IR spectrao support this conclusion and indicate that acyclic zig-zag chains are abundant in the liquid. 
Finally, the neutron scattering measurements!!!! of the radial correlation function of liquid DF yield H-F and F-F 
average neighbor distances consistent with a most probable local structure similar to that of the other forms of HF. 

Since all associated forms of HF share the (HF) 2 dimer as a common structural unit, any satisfactory potential model 
should reproduce the main features of the (HF)2 potential energy surface. Ab initio quantum mechanical calculations 
have contributed significawtly-,to clarify the interactions present in the isolated dimer (HF-W- by, determining accurately 
its equilibrium propertieaHic2l and by mapping out the entire potential energy surface .EJEjO Much of this work has 
been directed to,eb-taiB analytical models for the potential surface, usually by fitting the numerical results to a chosen 
functional form.OtiH23~c3 The an a ] yhi rALexpressions of these models, which have been developed mainly to study 
the rovibrational states of (HF) 2 ,EiJE3ul£II are usually very complex and therefore not suitable for simple MD or MC 
simulations. 

Fortunately, as already mentioned, there is evidence that the main features of the real molecular interactions 
may be app^ximately reproduced by much simpler empirical models with effective potentials. In fact, it has been 
long knownlij~E3 that the bent structure of the (HF) 2 dimer is dominated by the classical two-body electrostatic 
interactions between the permanent multipole moments of the monomer (mainly the dipolar and quadrupolar terms). 
The incorporation of many-body polarization effects gives rather small refinements on the predicted equilibrium 
angles,E2l a finding consistent with the experimental observation that the dipole moment of the dimer is only weakly 
enhanced relative to that of the isolated monomers 

The effects of the environment on the HF molecules do not seem to be large. However, a comparison of the structures 
in the sequence "monomer — > dimer — > hexamer — > liquid — > solid" (as shown by Tables [n] and ^ below) evidences 
a tendency in which a larger degree of association implies a small increase of the H-F bond distance, together with a 
much larger decrease of the F-F distancejjfrom 2.8 A in the gas phase dimer to 2.49 A in the solid). Unfortunately, as 
stressed by Rothlisberger and Parrinello,t3 the available three-site modelsLju refer to rigid molecules and, consequently, 
they cannot reproduce the relaxation of the interatomic distances in going from the gas phase dimer to the condensed 
phases. 

The previous observations suggest a relatively clear picture of the ingredients required for a satisfactory potential 
model for HF. First, a charge distribution is required which approximately reproduces the first few multipolar moments 
of the HF molecule. Second, the H-F bond cannot be considered as rigid and mtra-molecular interactions must be 
included to allow for the observed variations of the H-F and F-F distances in the various aggregation forms of HF. 
Finally, further atom-atom interactions must be introduced to model the remaining non-Coulombic intermolecular 
forces. Having all this in mind, we have tried to construct a new three-site model by simultaneously using data on 
the gas, liquid and solid phases. For this purpose: (a) the molecules are not rigid, but the H-F bond length can vary, 
and, (bl the parameters are fitted to theoretical and-jexperimental data, including the ab initio structure_qf the HF 
dimer ,Efl the room temperature density of liquid DFliJ and the experimental structure of solid DF at 4 K.EEl 

For the structure of the dimer we have used thc_pi initio (HF) 2 potential energy surface developed by Bunker, 
Jensen, Karpfen, Kofranek and Lishka (BJKKL)oSlj The BJKKL calculations, in which the potential energy 
of the (HF) 2 complex has been determined for over 1000 different configurations, represent the most complete and 
accurate scan of the energy surface to date. These quantum mechanical results are in excellent agreement with the 
experiments on gas phase (HF) 2 . 

The decision to add to the fit some data on solid and liquid HF (more precisely, DF) has been taken because of 
the partial failure of our preliminary models fitted only to the ab initio surface. These models reproduced well the 
zig-zag (HF)oo chains characteristic of the crystal, but gane totally wrong inter-chain distances and, furthermore, 
did not agree withr-the. experimental density of the liquidO This behavior was to be expected. In fact, the (HF) 2 
potential surfacctSEaEJ only accounts for the basic HF-HF interactions within a chain, and obviously exclude the weak 
long-range interactions responsible for the distances between different (HF)^ chains. Since the density of the liquid 
is affected by interactions between distant pairs of HF molecules in relative orientations which are not sampled in the 
solid, it is also understandable that only by fine-tuning the long range potential it has been possible to reproduce the 
experimental density of the liquid. Finally, it must be mentioned that the addition of data on solid and liquid HF to 
the fit gave only a small deterioration of the agreement with the ab initio data on the dimer. This fact confirms that 
solid and liquid data add information on regions of the potential surface that are not sampled by the dimer. 
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II. METHODS AND CALCULATIONS 



A. Potential model 

The potential model is represented by intra- and inter-molecular parts: 

^F ra (nJ F ) = D e {1 - exp[-a(r HF - r e )}} 2 , (1) 

V -inter . t/Cou! ■ T/non- Coul 
AB — V AB V AB ) \ Z ) 

^non-Coul = exp( _ BFF rpF ) _ CpF rp 6 ( 4 ) 

bjkklE fitted the intra- molecular part of their own ab initio surface with a Morse potential, eq. (1). Here r e 
represents the equilibrium H-F distance, D e the dissociation energy and a is an effective range parameter. Because 
of their simplicity and accuracy, the BJKKL functional form and parameter values are adopted in this paper. 

The Coulombic interactions between molecules A and B are modeled through three point charges for each HF 
monomer, two at the nuclear positions i?n and Rf, and the third at a position Rx along the H-F bond. In eq. (3) 
qi and qj are the fractional charges on the ith site of molecule A and jth site of molecule B, respectively; r.y is the 
distance between these sites. The motion of the site X is constrained so that it remains at the same relative position 
along the bond, 

Rx = (3R F + (1 - 0)Rh, (5) 

where (3 is an adjustable parameter between and 1. The charges are +q at both the H and F nuclei, and — 2q at 
the third site to preserve neutrality. This three-site charge model is related to that of Refs. [lj-[|. By allowing for 
changes in the H-F bond length, tlie-pcesent model effectively accounts for a part of the polarization effects. For solid 
and liquid DF the Ewald's methodcjEi has been used to ensure complete convergence of the Coulombic interactions 
(which have an infinite range). 

The remaining non-Coulombic part of the inter-molecular potential is represented in a simplified way, using only 
a Buckingham "exp— 6" atom-atom interaction between the fluorines, eq. (4). This term is meant to represent the 
interactions between the electronic clouds around two far away atoms. Since the hydrogens in HF are essentially 
bare nuclei, it makes good physical sense to avoid atom-atom interactions involving them. As a matter of fact, no 
improvement in the quality of the fit is found by adding similar H-H and H-F non-Coulombic interactions. 

A rather well defined hierarchy of interactions may be identified in the chosen potential model. The length of the 
HF monomer is solely determined by the intra-molecular potential. The structure of the dimer is also influenced 
by the position of the charge site X and by the F-F equilibrium distance, i.e. by the position of the minimum of 
the exp— 6 interaction between the fluorines. Finally, the charge and the remaining properties of the exp— 6 model, 
mainly the strength of the long range attractive term Cff r^p, an?ec t the structure and density of solid and liquid 
HF. The presence of this hierarchy implies that small changes in the charge and in the long range attraction can be 
compensated by the remaining parameters to maintain the correct monomer and dimer structures. 

B. Potential optimization 

The potential mcidel contains five adjustable parameters, j3, q, Aff, Bff and Cff, and three parameters fixed at 
the ab initio values,E-3 r e , a and D e . 

In a first series of attempts, we fitted the present model (as well as other preliminary ones) only to the ab initio 
potential energy of the dimer. The \ 2 deviation between the model and the ab initio surface was computed, for each 
given combination of parameter values, with the same relative weightsL-used by BJKKLo The \ 2 was minimized 
by searching the parameter space with Nelder-Mead simplex methodE3 The resulting potential was then tested 
by computing some properties of the other associated forms of HF. In particular, the liquid phase was studied by 
isothermal-isobaric MD simulations, as described in the next Subsection. As discussed in the introduction, these 
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preliminary models fitted only to the ab initio surface gave unsatisfactory results, so that it was decided to add more 
data to the fit. 

For this purpose, the equilibrium geometries at T = K of the HF dimer, hexamer and crystal were determined 
as a function of the potential parameters by minimizing, with the WMIN program,E3 the total potential energy with, 
respect to the structural parameters. The deviations of the calculated geometries from the ab initio dimer structured 
and from the experimental DF crystal structured at 4 K are then added to x 2 , with weights subjectively chosen to 
make the contributions from the surface, the dimer and the crystal roughly equal. 

Since the new set of parameters, although more satisfactory, still did not give the correct density of the liquid, 
it become necessary to tune the long range interactions. In a further set of minimization runs, a range of q and 
Cff values was searched, again with Nelder-Mead method. The three remaining parameters /?, App and -Bff where 
determined as a function of q and Cff by fitting the ab initio and crystal data. Each complete set of five parameters 
was then used in a short MD simulation to determine thc.-enuilibrium density of the liquid at 293 K and to add to 
the x 2 the deviation from the relevant experimental value. til This method, which involves two nested fit procedures, 
was found to be reasonably efficient and converged to a minimum in about fifty cycles. 

The main problem encountered in the fit was the noise in the computed liquid density due to insufficient MD 
equilibration with our computer time constraints. To reduce this noise, the MD equilibration was done in parallel 
with the potential optimization, by accepting after successful x 2 minimization cycles the final configuration of the 
MD run as the initial configuration for the next run. With this strategy the current MD configuration was always 
the one with the best potential parameters so far. Since the parameters change rather slowly, the simulated system 
tended to remain close to equilibrium. No structural data on the liquid, beside the density, has been included in the 
fit. This rather drastic choice avoids the repeated calculation of equilibrated radial correlation functions, which would 
have required even longer MD runs then those needed for equilibrated densities. 

As a technical detail, it must be noticed that no attempt was done to embed potential surface calculation, MD sim- 
ulation, dimer and crystal energy minimizations and the two nested Nelder-Mead procedures into a single monolithic 
program, which would have been unmanageably complex. Separate programs, calling each other as distinct processes 
at the operating system levelJ^J were used instead. The two Nelder-Mead procedures, in particular, were actually a 
single program invoking a second copy of itself. No special changes were required for the surface, MD and energy 
minimization programs. 

The optimal set of parameter values is shown in Tab. [j| and represents a compromise among the best results that 
can obtained separately for the dimer, the crystal and the liquid. As usual, no special physical significance should be 
attributed to the potential parameters. In fact, because of the possible compensation among different terms in the 
potential model, alternative slightly different sets of parameters might have been used. The model may be simply 
regarded as a tool to reproduce the observed data and predict new results. 

C. Molecular Dynamics 

The MD calculations employed 500 deuterium fluoride (LIFXmolecules, in a cube with periodic boundary conditions, 
and, using Andersen's isothermal-isobaric (NPT) method^ ~l£3 simulated a liquid sample in contact with a heat bath 
and subject to a hydrostatic pressure of I atm (« I0~ 4 GPa). The simulated liquid was obtained by melting and 
equilibrating at 293 K an initial crystalline configuration. The behavior of the system as a function of temperature 
was then determined by raising or lowering the bath temperature in steps of 10 K. Each temperature was maintained 
for at least 5 ps for equilibration and further 5 ps for analysis. The equations of motion were integrated using the 
velocity Verlet algorithm,EJEj~L3 with a time step of 0.25 fs. As previously described, each DF molecule carried three 
point charges, at D, F and X sites. Ciccotti, Ferrario and Ryckaert method for linear constraints^ has been used to 
maintain each massless X charge at a fixed fraction j3 of the DF bond. 

III. RESULTS 

The most important properties calculated with the present potential model for the mon omer, dimer, planar (HF)„ 
rings, hexamer, crystal and liquid forms of HF (or DF) a re co mpared in Tables ||, [Hi] and |v| with the available 
experimental and ab initio data. As described in section |II C| , the properties of the liquid have been determined 
through MD calculations, while the equilibrium geometries at K of the other forms of HF are found by minimizing 
the potential energy. 
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A. Monomer, dimer and cyclic polymers of HF 



The excellent results (Tab. ||) for the equilibrium bond length, dissociation energy antLspectroscopic parameters 
of the monomer, which all depend only on the parameters fixed to the BJKKL values,L3 indicate that the Morse 
model accurately reproduces the main features of the true intra-molecular potential. The spectroscopic parameters 
v e (harmonic frequency) and x e (anharmonicity constant) for the energy levels of a Morse oscillator, E n — hv e [{n + 
5) — x e (n + ^) 2 ],^re directly obtained from D e and a through v e = ay/ D e /fj, and x e — hv e /AD e , where /x is the 
reduced HF massO 

The multipole moments computed from our three-charges model (taking the molecular center of mass as origin) are 
very close to the experimental and ab initio moments of the monomer.acd The dipole and quadrupole values follow 
the trend of the ab initio results and therefore are slightly underestimated with respect to the experimental data. The 
octupole and hexadecupole moments are also in reasonable agreement with the ab initio calculations. This overall 
agreement is an indication of a good match between the electrostatic interactions in real HF and in the model. 

The computed minimum energy structure of the dimer (Tab. [n]) compares well with the experimental and ab initio 
data. □0QE3 The experimental F-F distance is excellently reproduced, as well as both the angles H'FF and FFH" of 
the bent equilibrium configuration. Moreover, the increased length of the HF molecule in going from the monomer to 
the dimer, and the slight length difference between the two HF intramolecular bondsji-j are well predicted. Since these 
length changes were the primary reason for allowing non-rigid molecules, such a behavior must be considered very 
satisfactory. Unfortunately, in spite of the excellent dimer geometry, the dimerization energy is clearly underestimated. 
This drawback was not completely unexpected. In fact, in real HF and in quantum mechanical models dimerization 
is accompanied by hydrogen bonding, which is not explicitly incorporated in the present classical treatment. 

As a further test of the potential, we have computed the equilibrium geometry of the planar (HF)„ rings, with 
n = 2, 3 ■ • • 8. As shown in Tab. [II, the structural par ameter computed for planar rings, with C n h symmetry, are in 
good agreement with the available ab initio results. Oc3 Though the model systematically underestimates the ab initio 
binding energies, it nevertheless reproduces correctly the relative stability of the different structures. The smallest 
ring, which is the cyclic dimer, has a binding energy of AE = —3.18 kcal/mole, and is thus substantially less stable 
than the bent dimer (Tab. ||). 

For cyclic (HF)3 the binding energy for each hydrogen bond, —AE/3, is slightly less that the binding energy of the 
bent dimer. The bond stabilization energy, —AE/n, increases for larger rings, up to the hexamer, and then decreases 
again. The particularly favorable stability of the (HF) n rings with n rj 6 is readily understood by noticing that for 
these rings the H'FF and FFH" angles (Tab. Ill) are close to those of the bent dimer (Tab. [Oj). Planar (HF) n rings 
must satisfy the geometric constraint FFH" — H'FF = a n , where a n = 180° — 360°/n is the inner angle of the n-sides 
regular polygon. The hexamer, for which a„ = 120°, can be obtained by joining essentially undeformed bent dimers, 
and is the most stable structure. 

The hexamer can be stabilized even further by allowing for non planar structures. We find that the stablest structure 
has a non-symmetriC|-^chair" shape, with average bond lengths and angles (Tab. |n|) which compare well with the 
available experimental!!! and which are almost identical to those found in the dimer. A "boat" structure slightly above 
in energy, at —4.96 kcal/mole, is also found, with lengths and angles close to those of the "chair" structure. 



B. Crystal and liquid DF 

Low temperature cpwstals of DF are orthorhombic, space group Cmc2i (C^ 2 ), with four molecules per unit cell 
on the o~ v plane sites £3 The minimum energy structure computed for the DF crystal (see Tab. IV) is close to the 



experimental structure at low temperature.Ej The discrepancies in the lengths of the cell axes partially compensate 
each other, yielding a density only slightly smaller than the experimental one. The increased H-F intra-molecular 
lengths with respect to both monomer and dimer, and the F-F distance smaller than in the dimer, are well reproduced. 

In Fig. [I] the densities predicted by the present model over the whole temperature range of liquid HF (which, at 
atmospheric pressure, goes from the freezing point at —83 °C up to the boiling point at 19.75 °C) are compared with 
the experimental data. The experimental densities show a nearly linear dependence on T. The data are from two 
sources, covering different temperature ranges.oH The slight vertical shift between the two data sets is due to the 
experimental uncertainties. It is very satisfactory to note that, although our potential model has been fitted only to 
a single density at 293 K, the straight line corresponding to its predictions is somewhat vertically shifted, but has 
essenti&UYnthe same slope as the experimental density. This slope is not reproduced by the HFC, HF-JV1 and HF-JV2 
models,B~Q since their density straight lines (obtained by joining the corresponding results, which, unfortunately, are 
available only at 203 and 273 K) intersect the experimental curve. 
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Another pleasant feature of the present model, also shown in Fig. [I], is that at 303 K and above the average density 
continued to decrease over the whole MD analysis period. At each T < 303 K, the density oscillated around an 
equilibrium value. In our opinion, this behavior indicates that the simulated liquid boils at some point in the range 
293 -T- 303 K, in good agreement with the experimental normal boiling point of HF (T& = 292.9 K). 

With regard to the internal energy U, our MD results exhibit a nearly linear dependence on T. However, their 
absolute values (Tab. IV) are underestimated with respect to the available experimental datau (a similar drawback is 
also present in the HF3 modellll). 

Fig. |^ reports the partial pair radial correlation functions gijij) computed for the liquid at 293 K, together with the 
real space function d{r) = 47rp m r [G(r) — 1]. Here p m is the molecular number density, and G{r) is a composite (or 
total) pair correlation function, obtained by adding the three partial pair correlation functions with the appropriate 



nuclear weig htsO and t 
scattering experiments 



^n by convoluting the sum with the same experimental resolution function used in the neutron 
The resulting theoretical prediction for d(r) is compared in Fig. |^, jrath the corresponding 
neutron diffraction data for liquid DF at 293 K, which are the only available real space data.EZI 

The first peak at r « 0.95 A is due to the intra-molecular H-F bond distance and here the agreement between 
simulation and experiment is excellent. The second and third peak of the experimental d(r) occur at r » 1.6 A 
and r ~ 2.55 A, and correspond to the hydrogen-bond thf inter-molecular distance and to the tff separation, 
respectively. Unfortunately, the model (with the present choice of parameters) fails to reproduce these peaks and 
the complex liquid structure at longer distances. The reason of such a shortcoming may be found by comparing the 
partial pair correlation functions gij(r) of the present model Fig. (Q) with the analogous ab initio MD results of 
Rothlisberger and Parrinello.Ea The position of the first inter- molecular peak of all our gij{r) is shifted toward larger 
r values (a similar trend occurs in the polarizable HF-JV2 modelo). In addition, the height of the hydrogen-bond 
peak of <7HF( r ) is nearly half the correct value, and the height of the first peak of ,9HH( r ) is also underestimated] 

It is to be recalled that the HF3 potential!! reproduces the three principal peaks of the experimental datatZl and 
gives the best performance for d(r) among the available models. This quite good agreement was obtained by modeling 
the hydrogen-bond interaction with a Morse termB To our knowledge, no other empirical model takes explicitly into 
account the hydrogen bond. 

Before concluding this Section, it may be useful to summarize the performances of the best available models for 
liquid HF. The HF3 potential gives, as already mentioned, quite good agreement with neutron diffraction structural 
data at the normal boiling point (293 K)JiH but systematically fails to reproduce thermodynamics^ the predicted 
internal energies are largely underestimated (i.e., their absolute values are too small) and the pressures in MD 
calculations at constant volume are several kilobars tooJiigh, indicating that the model system is less strongly bound 
than real HF. In comparison with HF3, the HFC modelcl yields a slightly better thermodynamics but a slightly worse 
structure for the liquid, with charges corresponding lo a dipole moment enhanced with respect to the monomer. 
Then, the predictions of HF-JV1 for the liquid phaseQ are not very different from the HFC ones, but HF-JV1 fails 
completely to reproduce the properties of the isolated dimer. Finally, whereas HF-JV1 works reasonably well only at 
room temperature, the polarizable model potential HF-JV2 represents a true improvement as regards the predicted 
density at low temperature, i.e. at 203 K. However, problems with the pair correlation functions are encountered also 
for HF-JV2. 



IV. CONCLUSIONS 



We have shown that a simple potential model suitable for MD or MC simulations can reproduce, quantitatively 
or semiquantitatively, many physical properties of the hydrogen fluoride over the whole set of its solid, liquid and 
gaseous associated forms. To cover such a wide range of environmental conditions, it has been necessary to consider a 
molecular model with variable bond length. The present investigation confirms the plausibility of assuming the dimer 
as the basic structural unit, but also stresses the need of fitting the potential to a set of experimental and theoretical 
data which includes information on the condensed phases. 

In the liquid, the correct trend of the density is very encouraging, but the underestimated energies and the problems 
with the radial distribution functions indicate that some physical effect is still misrepresented by the model. Unfortu- 
nately, it is not likely that adding data on the liquid structure could improve the results with the current type of model 
(without an hydrogen bond term). In fact, our experience with the fit shows that those few parameter sets which give 
better liquid structures were incompatible with the gas and crystal data. We believe that the essential shortcoming 
of the model is the neglect of any explicit representation for the hydrogen bond, which cannot be reduced to purely 
electrostatic multipole interactions. Such an explicit modeling of.±he hydrogen bond is also lacking in most other 
three-site potentials for HF, with the exception of the HF3 modelE Although the hydrogen bonding has a quantum 
mechanical origin, an approximate classical treatment is however possible and may have significant consequences, as 
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seen from the rather good structural results for the HF3 potential^ The inclusion of potential terms representing the 
hydrogen bond interactions appears therefore the next necessary step for more accurate HF models. More data on 
the liquid, including at least the radial correlation function, need to be incorporated in the fit. 

In conclusion, our model cannot be considered as a definitive solution of the problem, but can be seen as a significant 
step towards a really satisfactory potential for MD or MC simulations. It has the merit of pointing out the importance 
of a variable molecular length and of the hydrogen bond. Moreover, it shows that the strategy of a simultaneous fit 
to data covering all the associated forms of HF can be successful and should be considered as the appropriate way to 
fully accomplish the difficult task of finding a potential model for such a strongly associating system. 
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FIG. 1. Density of liquid DF as a function of temperature T. The center and half width of the error bars represent the 
average and standard deviation of the MD results, respectively, sampled over the last 5 ps at each bath temperature. The 
arrow at 303 K spans the range of densities sampled at this T, and indicates th*jt the density was decreasing through the 
whole analysisr-fi>p|iod. The filled circle represents the experimental DF density,E-J whereas the empty circles represent HF 
measurements,Ej'c3 multiplied by the DF/HF mass ratio to obtain consistent units. 



FIG. 2. Upper panel: total neutron pair radial correlation function d(r) for liquid DF (continuous line, simulation; dots, 
experimentO). Lower panel: gij(r) partial pair correlation functions from the simulation. The curves for <?hh(V) and <?FF(r) 
have been vertically displaced. 



TABLE I. Potential parameters. The meaning of the parameters is described in section [I A , The Morse parameters r e , a 
and D e have been fixed-to the values of the corresponding parameters (ci, C3 and kg 00 /8n 3/ ' 2 , respectively) of the BJKKL fit 
to the ab initio surfaceEj 

Monomer equilibrium length, r e (ao) 1.73727 

Monomer Morse parameter, a (a^ 1 ) 1.174 

Monomer dissociation energy, D e (Hartree) 0.22306 

Position of the X site, (3 = XF/HF 0.16245 

Charge, q, —2q, q on F, X, H, q (e units) 0.59456 

Buckingham parameter, Tff (kcal/mol) 167017 

Buckingham parameter, _Bfp (A -1 ) 4.148 

Buckingham parameter, CW (A 6 kcal/mol) 547.9 
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TABLE II. Properties of gas-phase HF monomer and polymers. The ab initio hexamer datallj are for a planar ring, whereas 
the model results are for the minimum energy "chair" structure. 



Monomer 



Model 



Experimental 



Ref. 



Ab initio 



Ref. 



HF equilibrium distance, r e (A) 0.9193 

dissociation energy, D e (kcal/mole) 140.0 

harmonic frequency, v e (cm -1 ) 4120.16 
anharmonicity constant, x e u e (cm -1 ) 86.69 

fundamental frequency, Auq^i (cm -1 ) 3946.78 

dipole moment, fi z (D) 1.772 

quadrupole moment, Q zz (D A) 2.122 

octupole moment, Q. zzz (D A 2 ) 1.894 

hexadecupole moment, & zzzz (D A 3 ) 1.658 

Dimer Model 



HF' distance (A) 

HF" distance (A) 

FF distance (A) 

H'FF angle (degrees) 

FFH" angle (degrees) 

dimerization energy AE (kcal/mole) 

Hexamer 



0.91680 
141.6 
4138.32 
89.88 
3961.42 
1.826 
2.36 



Experimental 



13 
13 

13 
13 
47 

41 
41 



Ref. 



0.9194 
141.2 
4135 
90.1 

1.7728 
2.3048 
1.7327 
1.87 

Ab initio 



Ref. 



0.9270 
0.9222 
2.6850 
7.52 
110.86 
-4.03 

Model 



2.72±0.03 
10±6 
117±6 
-4.56 

Experimental 



10 
10 

io 

42 



Ref. 



0.9236 
0.9220 
2.7919 
6.81 
114.45 
-4.32 

Ab initio 



Ref. 



HF distance (A) 

FF distance (A) 

H'FF angle (degrees) 

FFH" angle (degrees) 

binding energy AE (kcal/mole) 



0.9333 
2.6217 
6.15 
111.24 
-4.98 



0.973±0.009 
2.535±0.003 

104 
-7.20 



15 



15 
48 



0.948 
2.475 
2.4 
117.6 
-8.3 



TABLE III. Interatomic distances, angles, and bond stabilization energies of (HF) n planar rings with C„h symmetry. 









Model 








Ab initio 








n 


HF 


FF 


FFH" 


AE/n 


HF 


FF 


FFH" 


AE/n 


Ref. 




(A) 


(A) 


(degrees) 


(kcal/mol) 


(A) 


(A) 


(degrees) 


(kcal/mol) 






2 


0.9247 


2.6950 


49.34 


-1.59 


0.9223 


2.796 


54.23 


-3.30 


[ 




3 


0.9321 


2.6568 


81.37 


-3.79 


0.932 


2.616 


83.6 


-5.1 


[ 


iil 


4 


0.9335 


2.6246 


101.04 


-4.71 


0.943 


2.522 


101.6 


-7.2 


[ 




5 


0.9326 


2.6338 


114.41 


-4.90 














6 


0.9315 


2.6474 


123.90 


-4.92 


0.948 


2.475 


122.4 


-8.3 


[Ml 


7 


0.9308 


2.6599 


131.02 


-4.89 














8 


0.9302 


2.6688 


136.50 


-4.86 
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TABLE IV. Properties of solid and liquid DF 



Solid DF at 4.2K 



Model 



Experimental 



Ref. 



HF distance (A) 

FF distance (A) 

FFH" angle (degrees) 

unit cell axis a (A) 

unit cell axis 6 (A) 

unit cell axis c (A) 

density (g/cm 3 ) 

binding energy AE (kcal/mole) 



0.933 

2.662 

113.1 

3.29 

4.50 

5.31 

1.77 

-7.41 



0.95±0.02 

2.50±0.01 

116.6±1.0 

3.31 

4.26 

5.22 

1.89 

-10.7 



Liquid DF at 293K 



HF distance (A) 

FF distance (A) 

binding energy AE (kcal/mole) 

density (g/cm 3 ) 

boiling temperature Tt (K) 



Model 



Experimental 



Ref. 



0.930±0.022 
2.777±0.169 
-3.70 
1.02 ± 0.03 

< 303 



0.958±0.002 

2.56 

-6.93 

1.0106 

292.90 
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